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The role of Coulomb focusing in above-threshold ionization in an elliptically polarized mid-infrared strong 
laser field is investigated within a semiclassical model incorporating tunneling and Coulomb field effects. It 
is shown that Coulomb focusing up to moderate ellipticity values < 0.3) is dominated by multiple forward 
scattering of the ionized electron by the atomic core that creates a characteristic low-energy structure in the 
photoelectron spectrum and is responsible for the peculiar energy scaling of the ionization normalized yield 
along the major polarization axis. At higher ellipticities, the electron continuum dynamics is disturbed by the 
Coulomb field effect mostly at the exit of the ionization tunnel. Due to the latter, the normalized yield is found 
to be enhanced, with the enhancement factor being sharply pronounced at intermediate ellipticities. 

PACS numbers: 32.80.Rm; 32.80.Fb 



I. INTRODUCTION 

The rescattering of an ionized electron by the atomic core 
in a strong laser field plays a fundamental role in strong field 
physics d El, in particular, giving rise to above-threshold 
ionization, high-order harmonic generation and nonsequen- 
tial multiple ionization BHU- In a laser field of elliptical 
polarization, the ionized electron acquires a lateral drift mo- 
tion with respect to the major polarization axis (m.p.a.) which 
could seem to suppress the rescattering, with the consequence 
of extinguishing the related effects. However, recently it has 
been shown that, even in a laser field of elliptical polarization, 
the rescattering can play a significant role in nonsequential 
double-ionization l9lfT2l. Nonsequential double-ionization is 
due to back-scattering of the electron at small impact parame- 
ters. In contrast, the multiple forward- scattering of an ionized 
electron by the atomic core at large impact parameters induces 
Coulomb focusing (CF) fT3lfT5ll . Another contribution to CF 
comes from the initial Coulomb disturbance [the momentum 
transfer to the ionizing electron by a Coulomb field at the ini- 
tial part of the electron trajectory near the tunnel exit] [16]. 
Recent experiments |[T7lfT9l have shown that a characteristic 
spike-like low-energy structure (LES) arises in the energy dis- 
tribution of electrons emitted along the polarization direction 
of linearly polarized mid-infrared laser radiation. The mul- 
tiple forward scattering is shown to be responsible for LES 
l20H22l (for other aspects see (23 -26]). A question arises if 
the multiple forward scattering of the rescattering electron sur- 
vives in a field of elliptical polarization and, in general, how 
CF and LES are modified due to ellipticity of the field. 

Several experiments have been devoted to above-threshold 
ionization in an elliptically polarized laser field l27ti30l . in 
particular, the dodging phenomenon has been discovered [ 29 - 
l3TTl . It is expressed as a quick drop of the photoelectron nor- 
malized yield along m.p.a. with increasing ellipticity [the 
yield is normalized to the one in a linearly polarized field] that 
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FIG. 1 : (color online) The photoelectron normalized yield emitted 
along m.p.a. versus ellipticity and photoelectron energy: (a) Analyt- 
ical estimate of the yield in the PPT model via Eq. {7}; (b) Numer- 
ical simulations without taking into account Coulomb field effects; 
(c) Numerical simulations with Coulomb field effects. The laser in- 
tensity is Iq = 9 x 10 13 W/cm 2 and wavelength A = 2/im. The target 
atom is hydrogen. The curves with value 0.5 are for eye guide. 



is reversed when circular polarization is approached. The first 
indication of the Coulomb field effects for above-threshold 
ionization in an elliptically polarized field, manifested in the 
lack of the four-fold symmetry of the photoelectron angu- 
lar distribution, has been shown experimentally l28l and dis- 
cussed theoretically in 132H371 . 

In this paper, we investigate how the Coulomb field effects 
of the atomic core modify the dodging phenomenon in above- 
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threshold ionization in an elliptically polarized mid-infrared 
laser field. First, our attention is focused on the ellipticity 
and energy resolved photoelectron yield along m.p.a. up to 
moderate ellipticities. We show that due to CF, a remarkable 
energy dependence arises in the yield. What is more surpris- 
ing is the energy dependence of the yield for the low-energy 
domain [decreasing of the yield with increasing energy, in the 
energy interval of (0,4) eV in Fig. [TJ is reversed with respect 
to that for the high-energy domain [in the energy interval of 
(4,40) eV in Fig. [TJ. We demonstrate a direct relationship 
of this peculiar energy dependence of the yield with the LES 
appearance and with the specific features of multiple forward 
scattering of the ionized electron by the atomic core. Sec- 
ondly, we investigate the role of the Coulomb field effects at 
high ellipticities of the field. The CF is shown to enhance 
the photoelectron total yield along m.p.a. at high ellipticities. 
Moreover, the enhancement factor peaks at an intermediate 
value of ellipticity. 

We employ a semi-classical model incorporating tunneling 
and Coulomb field effects. Note that the CF effect and LES 
are conspicuous in mid-infrared laser fields when the Keldysh 
parameter is small, y = ^/l p /2U p C 1 0, where I p is the 
ionization potential and U p the ponderomotive energy. In this 
case the electron dynamics after tunneling is mainly classical 
and not disguised by quantum interference. 

The structure of the paper is as follows. In Sec. |Il|our theo- 
retical model is presented. Coulomb focusing up to moderate 



ellipticity values is considered in Sec. Ill In this section the 
appearance of LES and peculiar energy scaling of the normal- 



ized yield are shown and discussed. Sec. [IV] is devoted to the 
discussion of the Coulomb focusing effect at high ellipticities 
of the laser field. Sec. IVl concludes our discussion. 



II. THEORETICAL MODEL 

Our investigation is based on the classical-trajectory Monte 
Carlo method, with tunneling and the Coulomb field of the 
atomic core fully taken into account. The ionized electron 
wave packet is formed according to the Perelomov, Popov, 
Terent'ev (PPT) ionization rate [38 , 39 ] which is further prop- 
agated by the classical equations of motion: 



dp 

dt 



-E(f) + W c (r), 



(1) 



where Vc(r) — Z/r is the Coulomb potential of the atomic 
core [the target atom is hydrogen, Z = 1, and atomic units 
are used throughout the paper]. The laser field E(t) = 
(E x (t),E y (t),Q) is elliptically polarized: 



Ex(t) 

Ey(t) 



Eof(t) cos cot, 
—^Eof(t)sm cot, 



(2) 



with the ellipticity t; (\t;\ < 1). The electrons are born at the 
tunnel exit with the coordinates 



Xi - 



^cos/3, y;=x-sin/3, zi = 0, 



where j3 = arctan[— £ tan(jfy], cp t = cott is the ionization phase 
and x\ the initial position along the laser polarization direction 
(40). The initial momentum components are 



Pix = -mcosasin/3, 

Piy ~ Pi± cos a cos P •> 

Piz = A_Lsina, 



(4) 



where a is the angle between pt± and the axis y' [the y axis 
after rotation by an angle /3 around axis z] which is randomly 
distributed within the interval (0, lit). The transverse momen- 
tum pi± follows the corresponding PPT distribution |[38ll . The 
positions and momenta of electrons after interaction with the 
laser pulse are used to calculate the final asymptotic momenta 
PHI at the detector. Only electrons emitted along the x direc- 
tion (m.p.a.) within an opening angle Oo = ±2.5° are col- 
lected. The laser pulse profile is half-trapezoidal, constant 
for the first ten cycles and ramped off within the last three 
cycles. The number of propagated electrons is 10 6 and the 
convergence is checked via double increase of the electron 
number. The electrons are launched within the first half cycle 
(coti G [0, 7r]), since there are no multi-cycle interference ef- 
fects in the classical theory. The model has been confirmed to 
provide an adequate description for the strong field dynamics 
in the mid-infrared regime l20l . 



III. LOW-ENERGY STRUCTURE 

As we are concerned with the CF impact on the dodging 
phenomenon, let us first recall how this phenomenon arises 
when the Coulomb field effect is neglected l29l . The elec- 
trons born near the maximum of the elliptically polarized field 
tend to drift towards the minor axis of the polarization ellipse. 
Consequently, the electron should have a rather large trans- 
verse momentum \p y \ w |A y (<p^)| « %Eq/g) at the tunnel exit 
to counteract the drift and to reach the final state where the 
momentum points along m.p.a. (along x). Here A y ((pi) is the 
field vector-potential component at the ionization phase cp t . As 
a result, the ionization probability for electrons moving along 
m.p.a. is exponentially decreased with rising ellipticity (up to 



W, 



'^exp(-^E 2 /m 2 Al) 



(5) 



where 



/Hp is the width of the transverse momen- 



tum distribution 114211 . In describing the dodging phenomenon 
|29lk the photoelectron yield Y^(e) along m.p.a., normalized 
to the yield in a linearly polarized field, has been introduced 
and measured, where % and £ indicate the ellipticity and elec- 
tron energy dependence of the yield, respectively. Using 
the PPT probabilities for the photoelectron momentum dis- 
tribution in an elliptically polarized field d 3 W^/d 3 p 14211431 , 
where the CF effects for the electron dynamics in the contin- 
uum are neglected, the normalized yield can be derived: 



(3) 



r 5 (e)i 



La (d 3 W^/d 3 i>)d£l 

j AQ (d 3 w /d 3 v )da' 



(6) 
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Here the integration is carried out over the solid angle A£l 
along the m.p.a. with an opening angle 0n- At t; < y/E$ the 
yield reads 



75(e) = exp^-£ 



2 [ E l 

\ co 2 A 2 ± 



2V2e^E \ 
2 1 




where A| = 3£q/(27 p ) 3//2 &> 2 is the width of the longitudinal 
momentum distribution and I\ (Z) is a modified Bessel func- 
tion. As Eq. ([7]) shows, the yield exponentially decreases with 
rising ellipticity § and depends weakly on energy £, when the 
CF effects are neglected, see Fig. [T] (a, b). The cut across 
Fig. [I] at a fixed energy illustrates the first point, while the 
cut across Fig. [T] at a fixed ellipticity illustrates the second 
point [see also the dashed line in Fig. [5](c) below]. With the 
Coulomb field included, our numerical simulation shows an 
anomalous energy dependence of the yield for ellipticities up 
to ~ 0.3, see Fig. [T](c) and also Fig. [5](c) below: the yield 
decreases with increasing energy in the low-energy interval 
of (0,4) eV, while the energy dependence is reversed in the 
high-energy domain of (4,40) eV. 

To understand the role of the Coulomb field effects, we have 
calculated the photoelectron normalized yield, neglecting the 
momentum change induced by Coulomb field, either for the 
final longitudinal or for the final transverse momentum. This 
is accomplished by propagating an electron with certain ini- 
tial conditions twice: taking into account Vc(r) in Eq. {!]) or 
neglecting it. Afterwards, we replace the longitudinal (trans- 
verse) component of the exact final momentum with that de- 
rived from the calculation neglecting Vc(r). We noted that 
neglecting the Coulomb field for the longitudinal momentum 
has no significant effect on the yield, while the same for the 
transverse momentum, eliminates the peculiar energy depen- 
dence. As the CF is induced by the reduction of the transverse 
momentum of the ionized electron fT5lL we conclude that the 
CF is behind the anomalous energy dependence of the nor- 
malized yield. 

The anomalous energy dependence of the yield is closely 
connected with the appearance of LES in the energy spectrum 
of electrons emitted along m.p.a.. The latter is calculated nu- 
merically and shown in Figs. [2|4] We see that the LES [the 
peak in the energy interval (0, 4) eV] persists up to intermedi- 
ate ellipticities - 0.3 at 7 = 4.5 x 10 13 W/cm 2 , £ - 0.25 at 
7o = 9 x 10 13 W/cm 2 and £ - 0.23 at7 = 1.2 x 10 14 W/cm 2 ]. 
Figs. |2]|4]show also that at larger ellipticities, although there 
is no LES, the yield is significantly enhanced due to CF. The 
LES and CF due to soft forward scattering persist up to higher 
values of ellipticity at a given laser intensity than hard recolli- 
sions with backward scattering. The upper value of ellipticity 
for back- scattering (<^) can be estimated (9J [TO) by equating 
the lateral drift momentum %Eo/(Q to the width of the trans- 
verse momentum distribution A^: ^ ~ CO/ ' ^/Eq(2I p ) 1 ^ . For 
the parameters chosen in Fig. |5J ^ « 0.098, while the LES 
for the same parameters persists up to t; ^0.3. Nevertheless, 
the scaling with a laser intensity of the upper limit of ellip- 
ticity for LES is surprisingly well reproduced by the above 
estimation t; ~ 7 " 1 ^ 4 . 
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FIG. 2: (color online) Photoelectron spectra emitted along m.p.a. in 
an elliptically polarized field with different ellipticity £ , for the cases 
with CF taken into account (solid) and without CF (dashed). LES is 
observable up to t; « 0.25. The laser and atom parameters are the 
same as in Fig. [T] 
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FIG. 3: (color online) Same as Fig. [2] with a laser intensity of Io - 
4.5 x 10 13 W/cm 2 . LES is observable up to § w 0.3. 
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FIG. 4: (color online) Same as Fig. [2] with a laser intensity of Iq - 
1.2 x 10 14 W/cm 2 . LES is observable up to § w 0.23. 
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To explain the reason of the anomalous energy dependence 
of the normalized yield at moderate ellipticities, we proceed 
with an estimation of the yield taking into account CF. First, 
we express the momentum distribution of the ionized elec- 
trons at the detector d 3 W^/d 3 p via the distribution function 
over the transverse momentum (with respect to the field 
at the ionization moment) and the ionization phase (pi at the 
tunnel exit: 



d 3 W^ = d 3 W^(p ilVil _) 
d 3 p d(pid 2 p i± 



(8) 



where J? 



d(ffnP/±) • transformation Jacobian. The fi- 

d(P||>P±) 

nal momentum of the electron p is a function of the variables 
at the tunnel exit p = P£(<P*,P;_l) which depends on £ and 
is essentially disturbed by CF. Consequently, the dependence 
of the distribution function d 3 W^ (<P/(p) , p^ (p)) /d(pid 2 \) i± on 
the final momentum is also modified due to CF. Fortunately, 
this modification is negligible when the transverse momen- 
tum induced by CF 8p c ^ is rather small: Sp£ <C 



[derivation of this condition is given in Sec. [iVj see Eq. ( [26] ) 
below]. This is the case in the region of the anomalous en- 
ergy dependence of the yield, see e.g. Fig. [7](b,d). Thus, in 
this region the CF effect is mainly contained in the Jacobian 
. With this approximation, the normalized yield with CF 
effects Y^(e) can be expressed via the yield without CF Y% (e) 
given by Eq. |7j): Yg(e) « Y^(e)(J^ //£)(/ /^), where the 
index "C" indicates incorporation of the CF effects. The ra- 
tio of the Jacobians can be expressed by the ratio of the ini- 
tial momentum-space at the tunnel exit at a fixed momentum- 
space at the detector: 



f 
f 



d 2 V c 



(9) 



i±0 



Moreover, without CF effects J^/Jo^ d 2 p i± ^/d 2 p i±0 = 1, be- 
cause when CF is neglected, the final momentum- space at the 
detector equals that for the contributing electrons at the tunnel 
exit for any ellipticity. In conclusion, the yield with CF effects 
Y^(e) can be expressed via the yield without CF Y^ (e) and the 

ratio of the initial momentum- space volumes d 2 \) ( r ± ^/d- 
at a fixed final momentum- space: 



2 C 

Pllo 



Jf(£) 



(10) 
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FIG. 5: (color online) (a) The initial phase-space volume of con- 
tributing electrons d 2 p^ versus photoelectron energy at B, = and 

t; =0.15, numerical calculation (filled circles) and estimation (see 
text) (open circles); (b) The ratio of the initial momentum- space vol- 
umes d 2 p ( r ± ^/d 2 pf 1 _ , numerical (filled circles) and estimated (open 
circles); (c) The normalized yield, numerical (filled circles) and es- 
timated via Eq. |T0| ) (open circles) [is obtained by multiplying the 
numerical values of the volume ratio from Fig. |5|(b) by (e) of Eq. 
(jTJ)]. The dashed line shows (e) via Eq. jTj). The laser and atom 
parameters are the same as in Fig. [T] 



CF originates from the ratio of the momentum- space volumes 
of contributing electrons at the tunnel exit and is due to CF 
influence. 

Now, let us see why the ratio of the momentum- space vol- 
umes has a minimum. We can estimate the initial transverse 
momentum-space volume at the tunnel exit for the electrons 
which at the detector move along m.p.a., with a given energy 
interval as follows. In the case of linear polarization, the elec- 
tron gets a transverse momentum kick 5 due to CF during 
its motion from the tunnel exit to the detector, therefore, only 
the electrons which are distributed at the tunnel exit on a ring 
in a transverse momentum- space with a radius 8p c L will move 
along m.p.a. at the detector, see Fig. |6](a,b). The width of 
the ring is of the order of the final transverse momentum al- 
lowance 8pf± ~ pQ , with p = y/2e. Therefore, in the linear 
polarization case the initial momentum- space of contributing 
electrons is proportional to the transverse momentum change 
due to CF 



We have calculated numerically the ratio of the momentum- 



space volumes d 2 pf L( e /d 2 \)f ±0 , see Fig. pl(b), 

fTob [op( 



, and checked the 

validity of the estimation of Eq. (10) [open circles in Fig. [5] 
(c)] comparing it with the exact numerical calculations [filled 
circles in Fig. [5](c)]. One can see that in the energy interval of 
(2, 30) eV, which includes the region of peculiar behavior of 
yf (£), the above estimation describes the energy dependence 
of the yield appropriately. As Y% (e) without CF [dashed line 
in Fig. [5](c)] depends on energy monotonously, we can con- 
clude that the minimum at £ w 4 eV in the yield Y^(e) with 



<* 2 pf±o< 



' InpQodp^ 



(11) 



In an elliptically polarized laser field, the momentum- space 
distortion is more complex, see Fig. [6] (c,d). In this case, 
the transverse phase- space of contributing electrons at the tun- 
nel exit can be expressed via the components of the momen- 
tum change &Pyg induced by CF, see Fig. pne): ^ 2 P^ ~ 



2R(j)8p f± , with R = [(Sp^) 2 + (8p^) 2 ]/8p^ and sin</> 



5 



8 P%/ R > yielding 



£=0.15 



(^) 2 



8 p c 



sin 



(*/&) 2 - 



(12) 

In the most essential range of the photoelectron energy, see 
Fig. |7j 5/?^ <C 5/?^ which allows simplification of the esti- 
mate d 2 ^ ± ^/d 2 pf ±0 oc Sp^/Sp^Q. These estimations are in 
qualitative agreement with the exact calculations, as Figs. [5] 
(a,b) show. Therefore, the minimum in the ratio of the ini- 
tial phase-space volumes d 2 P^/d 2 P/io * s connected with 
the minimum of the transverse momentum change due to CF 
5 p c ^ . The latter for the cases of linear and elliptical polar- 
izations are calculated and shown in Fig. [7] using the method 
described in l22l . The partial contributions of high-order scat- 
tering events into the total momentum change are shown sep- 
arately. The total momentum change due to CF, in both cases, 
exhibits a pronounced minimum at the energy corresponding 
to the minimum of the ratio of the momentum- space volumes, 
see Fig. |5](b). This minimum marks the LES region and cor- 
responds to the threshold of the multiple forward scattering. 
The ionized electrons with final energies larger than 4 eV re- 
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FIG. 6: (color online) The initial [at the tunnel exit, marked by "i"] 
and the final [at the detector, marked by "f"] transverse momentum 
distributions of electrons with specified energies, (a,c) £ = 3 eV and 
(b,d) £ = 10 eV for the cases of linear (a,b) or elliptical polariza- 
tion = 0.15) (c,d), respectively The initial distribution shifted by 
A y (q>i) is indicated as "j 7 "; (e) To the estimation of the momentum- 
space volume, see text Eq. (12) . The laser and atom parameters are 
the same as in Fig. [T] 
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FIG. 7: (color online) The transverse momentum change due to CF 
versus photoelectron energy. The total value according to numeri- 
cal simulations (black line). The estimated contributions from high- 
order forward scattering (here up to 4th) and from initial CF (ICF) at 
the tunnel exit are shown separately. (a,c) £ = 0; (b,d) £ = 0.15. The 
laser and atom parameters are the same as in Fig. [T] 



scatter by the atomic core only once, while the electrons with 
smaller energies can re-scatter several times, inducing larger 
momentum change, see Fig. [7] It is remarkable that contribu- 
tion of the second forward scattering is more pronounced in 
the case of elliptical polarization than in the linear one, noted 
also in O. The high-order forward scatterings have an es- 
sential contribution to the transverse momentum change of 
the low energy photoelectrons. They are responsible for the 
creation of the minimum in the ratio of the momentum- space 
volumes and, consequently, for the creation of the minimum 
in the normalized yield at photoelectron low energies. An im- 
portant conclusion is that the multiple forward scatterings play 
a significant role in CF also in the elliptically polarized field. 
It is the reason why the LES persists up to moderate elliptici- 
ties. 



IV. COULOMB FOCUSING AT LARGE ELLIPTICITY 

At large ellipticities t; > \/Eq, the character of CF is essen- 
tially different. In this domain, the electron drift momentum 
is dominated by the vector-potential at the ionization moment, 
with a minor contribution from the initial transverse momen- 
tum [42]. The contribution of the high-order forward scat- 
tering in CF, as Fig. [SJa) shows, becomes monotonous in 
(Pi and perturbative. Moreover, the contribution of the for- 
ward scattering in CF becomes negligible relative to that of 
the initial Coulomb disturbance. The momentum transfer at 
the initial Coulomb disturbance is directed along the field at 
the ionization moment fT6lL i.e. it is perpendicular to the drift 
momentum, causing rotation of the initial momentum distri- 
bution around the axis perpendicular to the polarization plane, 
see also (5J. Therefore, the momentum- space volume will not 
be essentially modified by CF in this case. In fact, as one 



6 




Volume ratio 



Energy (eV) 

FIG. 8: (color online) (a) The transverse momentum change due to 
initial CF and the first forward scattering for large ellipticity values 
indicated in the inset, (b) The ratio of the momentum-space volumes 
a? 2 P/l<s;/tf? 2 P/lo versus photoelectron energy and ellipticity.The laser 
and atom parameters are the same as in Fig. [T] 



can see in Fig. [8](b) the ratio of the initial to final transverse 
momentum-space volumes is approximately one at large el- 
lipticities % > 0.25. 

The effect of rotation of the momentum distribution due to 
CF can be monitored by observing the enhancement of the 
photoelectron yield along m.p.a.. In other words, we investi- 
gate the modification of the dodging phenomenon due to CF. 
We calculate the total normalized yield along m.p.a., see Fig. 
[9] and its enhancement factor due to CF, see Fig. 10 We find 
that the enhancement factor is sharply peaked at an interme- 
diate ellipticity value £ m w 0.4, see Fig. 10 This value corre- 
sponds to the minimum of the yield along m.p.a., cf. with Fig. 
[9]E1. 

In the following, we give a simple analytical estimate of the 
CF impact on the tunneling probability that elucidates the ori- 
gin of the peak of the yield enhancement factor. The tunneling 
probability of an electron can be estimated with exponential 
accuracy 



Wi((phPi±) « exp (-2E a /3E((pi)-pj ± /A 2 A 



(13) 



where % = (Oti is the laser phase at the tunneling moment, pi± 
the initial transverse momentum with respect to the instanta- 
neous direction of the field, and E a = (2I p ) 3 / 2 the atomic field. 
The electron momentum at the detector is 



p = p ^ •-A(<jP ^ •) + <5p , 



c 



(14) 



where p ; is the electron initial momentum at the tunnel exit, 
A(<pj) the vector-potential at the ionization moment and <5p c 
the momentum change due to CF [in the case of large elliptic- 
ities, it is due to initial CF]. 

In the limit £ — » 1 , the ionization of electrons which at the 
detector move along m.p.a. takes place mostly near the laser 



phase when the vector potential points along m.p.a., see Fig. 
11 (a). For the field given by Eq. ^ the vector-potential is 



M<Pd = ~ A o sin (pi, A y {(pi) = -£A cos(ft, (15) 

with Ao = Eq/cq, and the ionization phase is close to <j9; w 71 / 2. 
Introducing ft • = 71 / 2 — q)j ; <C 1 and expanding the field over a 
small parameter ft up to the second order, we obtain for the 
total field 



Eh 



^£o[l + 0, 2 (l-£ 2 )/2^ 



(16) 



The ionization phase is determined from vanishing of the final 
transverse momentum with respect to m.p.a., namely, p y = 
Pi y — A y ((pi) + 5 p^ « 0. The probability is maximal at pi y « 
0. Taking into account also that the longitudinal component 
(with respect to the laser field) of the momentum change due 
to initial CF is [16 ] 



8pC^KE((pi)/(2I p fl\ 



(17) 



we find that 8p^ k^Eq/ (2I p ) 3 / 2 and ft « 7i(o/(2I p ) 3 ^ 2 . As 
a result, the ionization probability in the limit § — ^ 1 is 



■f|g_>i °cexp< 



2E a 7T 2 G) 2 (l-g 2 ) ] 

'3^o 3E ^(2I p )V^}' 



(18) 



The last term in Eq. ([18]) describes influence of the initial 
CF on the ionization probability and leads to the probabil- 
ity enhancement. This is because the electrons which move 
along m.p.a., are tunneled out at a larger ft , consequently, at a 
stronger laser field, see Eq. ( [T6] ), when the CF effect is taken 
into account. 

In the opposite limit £ — » 0, ionization of the electrons 
which move along m.p.a. takes place mostly at (pi w 0. In 
fact, the ionization probability at the laser phase (pi « 71 j 2 is 



O 

> 

LL 
O 




1E-14 



Ellipticity £ 



FIG. 9: (color online) The photoelectron normalized yield along 
m.p.a. integrated over the total energy range for the cases with CF 
(Y^ F , black, filled circles) and without CF (Y^ CF , red, open circles). 
The laser and atom parameters are the same as in Fig. [T] 
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(b) 



- 2(im, 9x10 ,3 W/cm 2 
I — 2|im, 1.2x10 l4 W/cm 
/< 3nm,9x10 ,3 W/cm 2 

j-'A 



0.0 0.2 0.4 0.6 0.8 1.0 0.0 0.2 0.4 0.6 0.8 

Ellipticity C Ellipticity C 




(e) 



3nm, 
9x10 13 W/cm 2 



1.2 3 

1.1 CL 



0.0 0.2 0.4 0.6 0.8 1.0 

Ellipticity £ 



0.0 0.2 0.4 0.6 0.8 1.0 

Ellipticity C 



0.0 0.2 0.4 0.6 0.8 

Ellipticity C 



FIG. 10: (color online) The ratio Y CF /Y NCF for different laser in- 
tensities (To) and wavelengths (A): (a) numerical simulations, (b) 
analytical estimations, according to Eq. {251), (c) Iq = 9 x 10 13 
W/cm 2 , X = 2 jUm; (d) 7 = 1.2 x 10 14 W/cIrr 2 , A = 2 and (e) 
7q — 9 x 10 13 W/cm 2 , A = 3 jim\ numerical simulation (solid), ana- 
lytical estimation (dashed). 



negligible, because of weakness of the field at this phase (E ~ 
%Eq). Meanwhile, at (pi w 0, the field is large (E ~ Eq), the 
laser induced lateral drift momentum is small [as A y ~ £ Eq/co 
with < 1] and can be compensated by the initial transverse 
momentum, see Fig. 



11 (b). The total field expanded near 



(pi w reads 



E((pi)=E [l-(p?(l-t; 2 )/2}. 



(19) 



In this case, the initial transverse momentum p t _\_ is required 
to have a nonzero value such that pi y = pi± cos (pi would com- 
pensate the drift velocity and result in vanishing of the final 
lateral momentum. The momentum change due to the initial 
CF is SpC « -K^E (pi/(2I p ) 3 / 2 . Vanishing of the final lateral 
momentum of the electron p y & 0, leads to 



Pi± cos (pi 



w sln< * 



CO 



COS (pi. 



(20) 



From the latter pi± is derived employing an expansion over 
(pi- 



7TCO(Pi 



co l(2i p y/ 2 



1 



(21) 



Here, the first term between the brackets is due to the initial 
CF. The ionization phase (ft is determined by p x , which can 



be deduced from the x-projection of Eq. ( 14 ) 



<Pi 



Eo 

CO 



2pi±E 



(Up) 2 



(22) 



The latter simplifies for p x ^> tcEq [the electron energy larger 
than 0.3 eV for the chosen paremeters], yielding (pi « p x co/Eo. 



Inserting the obtained values for (ft and pi± into Eq. ( 13 ), one 
gets 



Wi(#c)|§ 



^0 



exp(-2E a /3Eo) 



x exp 



A 2 , 



A 2 



Eo 
(0 



1Z(OPx 



Integrating Eq. ( [23] ) over p x and keeping leading terms in the 
exponent (up to q^) over expansion by a small parameter t; , 
we derive the ionization probability in the limit £ — >• with 
exponential accuracy 



H>i|§_>oocexp 



(_2Ea 
I 3^o 



-2/72 
? ^0 



0) 



2 A 2 ' coi(2I p y/2j' 



(24) 



The last term in the exponent in Eq. ( [24] ) stems from the initial 
CF. It enhances the ionization probability because the initial 
transverse momentum of the electron, moving finally along 
m.p.a., is decreased when CF is accounted for, see Eq. pT 



Using the derived asymptotic expressions for the ionization 
probability given by Eqs. dT8l) and p4|), we interpolate the 



(25) 



ionization probability for arbitrary £ : 

■/(§) ~max{w;|^ 05 w*l<^i}> 



and from the latter estimate the ratio of the yield with to the 
yield without CF, see Fig. [TUJ The estimated ratio reproduces 
correctly the position of the peak. The CF effect is described 
by the last terms in the exponent in Eqs. ( [T8] ) and ( [24] ). These 
terms increase as one moves away from the asymptotic values 
of £ . Consequently, enhancement of the ionization probability 
due to CF is largest at an intermediate ellipticity |45|. The 
peak of the yield enhancement shifts to lower ellipticities with 
increasing intensity or wavelength. It is correctly predicted by 
our estimation via Eqs. ([T8]) and p4|), see Fig. 



[10] However, 

our qualitative estimates are carried out only with exponential 
accuracy, neglecting the pref actors. For this reason we cannot 
predict correctly the absolute value of the yield ratio. 

Thus, the normalized yield at large ellipticities is enhanced 
due to the Coulomb field effect [initial CF]. The reason for 
the enhancement is that the electrons which drift along m.p.a. 




FIG. 1 1 : (color online) The geometry of the ionization near the phase 
of the laser field when the ionization probability is maximal: (a) £ — > 
1, (b) § — >• 0. x-axis is along m.p.a. 
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have larger tunneling probability when CF is taken into ac- 
count. This is due to the fact that either [in the case of t; — » 1] 
the electrons are tunneled out at larger laser field values or [in 
the case of t; — » 0] the electrons are tunneled out with smaller 
initial transverse momenta when the initial CF is taken into 
account. At intermediate ellipticities both of the enhancement 
mechanisms can contribute which induces a sharp peak in the 
enhancement factor. 

Eq. ( [23) of this section can be used to deduce the condi- 
tion for neglecting the CF effects in the distribution function 
d 3 W^((pi(p),p i± (p))/d(pid 2 p i± at small The CF modifi- 
cation is given by the term proportional to 7lCQp x / (21 p ) 3 / 2 in 
the exponent of Eq. ( [23] ) which equals the momentum change 
due to CF 8p c ^. The modification is negligible if 8p c ^ <C 
and Sp^Eo/co <C Aj_. The second of these conditions is the 
strongest and reads 



(0 



For the parameters used, this yields 8p c ^ <C 0.14. 



(26) 



scattering of the ionized electron by the atomic core has dom- 
inated contribution in Coulomb focusing up to moderate ellip- 
ticity values. The multiple forward scattering causes squeez- 
ing of the transverse momentum- space volume, which is the 
main factor influencing the normalized yield at moderate el- 
lipticities, described by Eq. ( fTO) . It is responsible for the pe- 
culiar energy scaling of the ionization normalized yield along 
the major polarization axis and for the creation of a character- 
istic low-energy structure in photoelectron spectrum. 

At large ellipticities, the main CF effect is due to initial 
Coulomb disturbance at the exit of the ionization tunnel. The 
initial Coulomb disturbance, as our estimates in Eqs. ([18) 



and ( [24] ) show, enhances the ionization yield. This is because 
the electrons are tunneled out at larger laser fields or with 
smaller initial transverse momentum when the initial CF is 
taken into account for the electron drifting along m.p.a. The 
enhancement factor is shown to be sharply pronounced at 
intermediate ellipticities when both of the above mentioned 
enhancement mechanisms contribute. In this region of 
ellipticity, the yield is enhanced by an order of magnitude due 
to CF. 



V. CONCLUSION 

We have investigated the role of Coulomb focusing in 
above-threshold ionization in a mid-infrared laser field of el- 
liptical polarization. We have shown that multiple forward 
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and 



volumes d 2 p^ / d 2 pf^ discussed in Sec. I 
Fig.|8](b), while the estimates of Eqs. |T8] > and ( [24] ) of tETs sec 
tion deal with the CF influence on the tunneling probability at 
the tunnel exit and do not include the ratio of the momentum- 
space volumes. 



